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Abstract 

We use the cluster variation method (CVM) to investigate the phase structure 
of the 3d gonihedric Ising actions defined by Savvidy and Wegner. The geometrical 
spin cluster boundaries in these systems serve as models for the string worldsheets 
of the gonihedric string embedded in Z^. The models are interesting from the 
statistical mechanical point of view because they have a vanishing bare surface 
tension. As a result the action depends only on the angles of the discrete surface 
and not on the area, which is the antithesis of the standard 3d Ising model. 

The results obtained with the CVM are in good agreement with Monte Carlo 
simulations for the critical temperatures and the order of the transition as the self- 
avoidance coupling K is varied. The value of the magnetization critical exponent 
(3 = 0.062 it 0.003, calculated with the cluster variation-Pade approximant method, 
is also close to the simulation results. 



1 The Model 



In this paper we discuss the use of the cluster variation method (CVM) in mapping out 
the phase diagram of the gonihedric 3d Ising model. The gonihedric 3d Ising model is 
a generalization of the usual 3d Ising model where planar Peierls boundaries between + 
and — spins can be created at zero energy cost. It has been introduced in in relations 
with string theory. The CVM method in other contexts has shown itself to be an accurate 
and economical way of describing the phase diagram both for first order and continuous 
transitions. As both the model and the method may be unfamiliar, we outline both in 
turn before going on to describe our results. 

The genesis of the model is in a novel discretized random surface theory, the so-called 
gonihedric string 

s = W\x^-x,\e{a.^), (1) 

where the sum is over the edges of some triangulated surface, 9(aij) = \n — aij\'', ( is some 
exponent, and aij is the dihedral angle between neighbouring triangles with common link 
{ij). This definition of the action was inspired by the geometrical notion of the linear size 
of a surface, as originally defined by Steiner 0. 

In eq. (|l|) the surface itself is discretized, rather than the space in which it is embedded. 
An alternative approach to discretizing the linear size would be to restrict the allowed 
surfaces to the plaquettes of a (hyper)cubic lattice, which corresponds to also discretizing 
the target space. Savvidy and Wegner |^, ^ ^ did this and rewrote the resulting model 
as an equivalent generalized Ising model using the geometrical spin cluster boundaries 
to define the surfaces. The energy of a surface on a cubic lattice is then given hj E = 
n2 + 4:K,n4, where n2 is the number of links where two plaquettes meet at a right angle, 
77-4 is the number of links where four plaquettes meet at right angles, and k is a free 
parameter which determines the relative weight of a self-intersection of the surface. In 
the limit k —>■ oo the surfaces are strongly self-avoiding, whereas the opposite limit k ^ 
is that of "phantom" surfaces that pass through themselves without any energy penalty. 
It is worth emphasizing that the energy is very different from that of the standard 3d 
Ising model with nearest neighbour interactions where the surfaces are weighted entirely 
by their areas and not at all by their embeddings. 

On a cubic lattice the generalized gonihedric Ising hamiltonian which reproduces the 
energy E = n2 + ^nrii contains nearest neighbour next to nearest neighbour 

{{{hj))) and round a plaquette {[i,j,k,l]) terms 

-i7 = 2K^(Ti(Tj - ^ + (^i(^j(^k(Ti. (2) 

(ij) [id,k,l] 

Such generalized Ising actions and their equivalent surface formulations have quite com- 
plicated phase structures for generic choices of the couplings 0, ||, 0. The particular 
ratio of couplings in eq. however, is special and introduces a novel symmetry into the 
model, related to a zero-temperature high degeneracy point where it is possible to flip 
any plane of spins at zero energy cost. 
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Various approaches have been used to investigate these models, including a zero- 
temperature analysis, mean-field theory and Monte Carlo simulations. We briefly outline 
the results thus obtained for comparison with our CVM calculation in this paper. 

In an analysis of the zero temperature ground states of the model we write the full 
lattice hamiltonian as a sum over individual cube hamiltonians he and observe that if 
the lattice can be tiled by a cube conflguration minimizing the individual he then the 
ground state energy density is cq = min he [0. This approach conflrms that a layered 
ground state with parallel layers of flipped spins perpendicular to one of the lattice axes 
and arbitrary interlayer spacing is degenerate with the ferromagnetic ground state for all 
K. In addition, at k = extra ground states appear where all the spins of one of the 2 
sublattices of the cubic lattice are flipped 0, |TD . 



In the mean fleld approximation the spins are replaced by average site magnetiza- 
tions. The calculation of the mean fleld free energy can be performed by considering 
independent magnetizations for the 8 sites of an elementary cube. Therefore, as in the 
zero temperature approach, the energy can be still decomposed into a sum of individ- 
ual cube terms. Numerical iteration of the resulting eight coupled mean-fleld equations 
shows a single transition from a paramagnetic high temperature state to a layered, or the 
equivalent ferromagnetic, low temperature state The inverse critical temperature (3c 



determined in this fashion decreases quite sharply with k. 

The flip symmetry of the model poses something of a problem when carrying out Monte 
Carlo simulations. A simple ferromagnetic order parameter such as the magnetization 



^= Tl^E^O (3) 




will be zero in general, because of the layered nature of the ground state. Staggered mag- 
netizations also fail as order parameters because the interlayer spacing can be arbitrary. 
In [|rU| boundary conditions have been suitably chosen in order to pick out the ferromag- 



netic ground state allowing the use of standard (unstaggered) magnetization to extract 
magnetic critical exponents. Monte Carlo simulations with such boundary conditions 
||T0| for different n values on lattices of various sizes allowed a flnite size scaling analysis 
to be carried out in order to extract estimates for some of the critical exponents. For 
K = 1 this gave v = 1.2(1) from the ratios of slopes of Binder's magnetization cumulants, 
7/1/ = 1.79(4) from the FSS of the susceptibility x, and (a — l)/z/ = —1.3(2) from the FSS 
of the energy with two different sorts of flxed boundary conditions. All these exponents, 
rather remarkably given that the model is deflned in three dimensions, are close to the 
Onsager values of the two-dimensional Ising model with nearest neighbour interactions, 
as was the critical temperature Pc = 0.44. 

Simulations of other n > 1 values gave very similar results. However, the k = model 
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H = ^ (^i(^j(^k(ri. (4) 

[i,j,k,l] 

appeared to be a special case, displaying a flrst order transition. The transition stayed 
flrst order at k = 0.1 but softened rapidly as k increased, so the crossover to the second 
order behaviour seen at k = 1 was quite sharp. 
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2 The Cluster Variation Method 



The cluster variation method, or CVM for short, is based on a truncation of the clus- 
ter (cumulant) expansion of the free energy density functional on which the variational 



formulation of statistical mechanics is based ||12| , p!3[] . Unlike mean field theory it gener- 
ally locates rather accurately the boundaries between different phases in complex phase 
diagrams and, using the recently proposed cluster variation-Pade approximant method 



I^ , |T3|, one can extract non-classical, precise estimates of the critical exponents. 



For a generic Ising-like model described by a hamiltonian if on a lattice A the exact 
free energy can in principle be obtained by minimizing the functional 

F[pj,] = TT(^pj,H+^p^\npj,y (5) 

of the trial density matrix pA, subject to Tr(pA) = 1. 

Assuming that the hamiltonian can be written as a sum of cluster contributions 

H=Y.ha, (6) 
where F is a collection of clusters that suffices to enumerate the interactions, an approxi- 



mate free energy functional can be written in the form |13 



F = ^ Tr (p, + i ^ a,Tr(p„lnp,) (7) 

where P is a suitable set of clusters (which must contain F as a subset), the largest of 
which (called maximal clusters) have to reflect in some way the symmetry of the lattice, 
and the are numerical coefficients determined by 

^ = 1, Va G P. (8) 

The constraints are now 

Tr(p,) = 1, a e P 

pa = TTf3\a{pf3), acp (9) 

and the latter can actually be used as a definition of the density matrices of the subclusters 
of the maximal clusters. 

An important feature of the CVM is that the local minima of the approximate free 
energy can be easily found by means of a simple iterative procedure called natural iteration 
method ||T7| , |T^ , which has the property that, for any given initial set of density matrices, 
the iteration always converges to a local minimum of the free energy. 

In the present work, we have used the cube approximation of the CVM, that is the 
maximal clusters of the set P are the elementary cubic cells of our simple cubic lattice; 



for this approximation the free energy density functional has the form (see |jT2| for an 
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application to the simple nearest neighbour Ising model) 



f[ps] = TTipsHs) + ^ 



TtC^Ps) - ^ Tr£(p4_piaquette 
plaquettes 



-\ Tr/:(p2,edge) Tri:(pi,site) 

edges sites 



(10) 



where is the contribution of a single cube to the hamiltonian (when splitting the 
total hamiltonian H into single cube contributions one has to keep in mind that nearest 
neighbour interactions are shared by four cubes and then will get a coefficient 1/4 in Hg, 
and similarly next-nearest neighbour and plaquette interactions will get a coefficient 1/2), 
jC{x) = a:lna:, pa with a = 8 (4, 2, 1) denotes the cube (respectively plaquette, edge, site) 
density matrix, and the sums in the entropy part are over all plaquettes (edges, sites) of 
a single cube. Notice that we have not assumed any a priori symmetry property for our 
density matrices. 

The cluster variation method can be viewed as a generalized mean field theory, and 
hence it is clear that it can give only classical predictions for the critical exponents. In 
order to overcome this difficulty, one can use the recently proposed cluster variation-Pade 
approximant method (CVPAM) |]T^ |T^, which has proven to be a rather accurate 
technique, although not very demanding in terms of computer time. The basic idea of 
the CVPAM is that, since the CVM gives, for Ising-like models, very accurate results 
at low and high temperatures (i.e. far enough from the critical point), one can try to 
extrapolate this results in order to study the critical behaviour. In order to determine the 
critical exponent of the order parameter m, for example, one calculates m{(3) with the 
CVM up to a temperature at which the error can be estimated to be very small (typically 
of order 10~^), and then constructs, by a simple interpolation, Pade approximants for 
the logarithmic derivative of m{j3): the pole and the corresponding residue of each Pade 
approximant are then estimates for the critical temperature and for the critical exponent 
respectively. 



3 The Results 

In this section we describe our results for the phase diagram of the model (^ at different 
values of the parameter k. When the temperature is lowered, the model (|^) undergoes 
a phase transition towards a low temperature ordered phase: this transition has been 
investigated by means of the CVM in the cube approximation described in the preceding 
section. 

First of all we describe our results at k = 1. We have solved numerically our approxi- 
mate variational principle for j3 ranging in the interval < /? < 0.5 and we have calculated 
the values fp{/3), friP) and /l(/3) of the local minima of the free energy corresponding 
respectively to the paramagnetic, ferromagnetic and layered phases. The three functions 
/p(/3), fF{f3) and hiP) are plotted in Fig. |l|. 

When f3 < 0.404 the unique local minimum is that corresponding to the paramagnetic 
phase; in the interval 0.404 < (3 < 0.427 we find both the paramagnetic and layered 
minima of the free energy and fp{(3) < fiiP) in the whole interval. At low temperature 
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when (3 > 0.427 the ferromagnetic local minimum appears and it happens to be the global 
minimum of the free energy. Therefore the CVM predicts that the zero temperature 
degeneracy between phases with all possible sequences of "+" and "— " planes is broken 
at finite temperature. 

The transition at (3c = 0.427 from the paramagnetic to the low temperature ferromag- 
netic phase is critical. Indeed, when the inverse temperature j3 is lowered below Pc, the 
ferromagnetic local minimum of the free energy disappears. 

We have studied the model (|^) for other values of the parameter k ranging in the 
interval < k < 10. For each k > we have found that at low temperature the model is 
in the ferromagnetic phase and we have calculated the inverse transition temperature jSc] 
in Fig. El we have plotted (3c function of k. 

At sufficiently low values of n, that is k < = 0.87±0.01, the nature of the transition 
changes over to a first order behaviour which is strengthened as k is lowered. 

In Table |1] our results for (3c are compared with Monte Carlo and mean field approxi- 



mation results obtained in 101; Monte Carlo and CVM predictions are in good agreement. 



Table 1 



0.25 0.5 
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(3c CVM 0.550 0.464 0.443 0.427 0.421 0.420 0.420 

(3c MC0' P 0.505 - 0.44 0.44 0.44 0.44 0.44 



(3c MF0 0.325 0.31 0.278 0.167 0.09 0.0335 0.02 



Table 1: Our results for the inverse transition temperature (3c at different values of the 
parameter k are listed. These CVM results are compared with Mean Field and Monte 
Carlo previous results. 



At K = the first order transition is at (3c = 0.550. Moreover at low temperature the 
layered phases, the ferromagnetic phase and their antiferromagnetic versions obtained by 
fiipping the site-magnetizations of one of the 2 sublattices of the cubic lattice coexist with 
the same free energy. In Fig. ^ the free energy of the paramagnetic phase (dash-dotted 
line) and of the coexisting ordered phases (solid line) are depicted. 
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Therefore, at n = the zero temperature symmetry of the model is not violated, and 
the CVM correctly respects the ferro-antiferromagnetic exact symmetry of the partition 
function. 

Finally, for some values of k such that the transition is critical, that is k = 1, 2, 5, 10, we 
have evaluated the magnetization critical exponent (3 by means of the cluster variation- 

H, |l6l, obtaining 0.059 < /3 < 0.062 at k 
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1 and 



Fade approximant method 
P ~ 0.065 at K = 2, 5, 10. These results clearly suggest that the exponent is independent 
of K and a prudential estimate is (3 = 0.062 ± 0.003, which has to be compared with the 
conjectured Onsager value 1/8 on one side, but also with the estimate, based on finite 
size scaling of Monte Carlo results |jTO|, = 0.04(1): given the Monte Carlo estimate 
u = 1.2(1) one can say that the /? values predicted by CVM and simulations are in rather 
remarkable agreement. 



4 Conclusions 

In this work we have applied the cube approximation of the Cluster Variation Method 
to find the phase diagram of the gonihedric 3d Ising model defined in eq. (^. Moreover 
the low-temperature CVM results for the order parameter have been used to evaluate the 
exponent /3 at k = 1, 2, 5, 10 via the cluster variation-Fade approximant method. 

We summarize here our main results. The CVM approximation gives values of the 
inverse critical temperature in quite good agreement with those predicted by Monte Carlo 
simulations. The transition remains critical for k, > Utr = 0.87±0.01, where it becomes of 
first order. Our evaluation of the critical exponent j3 = 0.062 ±0.003 of the magnetization 
is also close to the simulations results of |jlO| . 



The new result of this paper is that the CVM predicts at finite temperatures a violation 
of the symmetry of the hamiltonian (^. We find that in the ordered region of the model 
the ferromagnetic phase is always stable with respect to the lamellar phase. This suggests 
to study the model (Ef) in a parameter space larger than the one used in this paper. The 



knowledge of the global topology of the phase diagram in an enlarged parameter space 
could be useful to answer the main question set by the Monte Carlo results and confirmed 
by the results of this paper about the nature of the critical transition of the model (0) at 
sufficiently high values of k. 
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Figure 1: Solid, dash-dotted and dashed hnes represent respectively the free energies 
/F(/9),/p(/3)and^(/3)at« = l. 
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Figure 2: The inverse transition temperature is plotted as a function of the parameter 
K. Dashed and solid lines represent respectively first and second order inverse transition 
temperatures. 
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Figure 3: The dash-dotted and the sohd hnes represent respectively the free energy of the 
paramagnetic phase and of the coexisting low temperature phases at k — 0. 
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